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ABSTRACT 

We investigate the relation between star formation rates and local galaxy envi- 
ronment for a stellar mass selected galaxy sample in the redshift range 1.5 < z < 3. 
We use near-infra-red imaging from an extremely deep Hubble Space Telescope sur- 
vey, the GOODS-NICMOS Survey (GNS) to measure local galaxy densities based 
on the nearest neighbour approach, while star-formation rates are estimated from 
rest-frame UV-fluxes. Due to our imaging depth we can examine galaxies down to a 
colour-independent stellar mass completeness limit of log M* = 9.5 M Q at z ~ 3. 
We find a strong dependence of star formation activity on galaxy stellar mass over 
the whole redshift range, which does not depend on local environment. The average 
star formation rates are largely independent of local environment apart from in the 
highest relative over-densities. Galaxies in over-densities of a factor of > 5 have on 
average lower star formation rates by a factor of 2 — 3, but only up to redshifts of 
z ~ 2. We do not see any evidence for AGN activity influencing these relations. We 
also investigate the influence of the very local environment on star-formation activity 
by counting neighbours within 30 kpc radius. This shows that galaxies with two or 
more close neighbours have on average significantly lower star formation rates as well 
as lower specific star formation rates up to z ~ 2.5. We suggest that this might be due 
to star formation quenching induced by galaxy merging processes. 
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1 INTRODUCTION 

Observational studies have shown that up to half of the cur- 
rently existing stellar mas s was already in p l ace by z ~ 1 
terinchmann fe Ellis! |2000|; iDrorv et alj 12004: iBundv et all 



Mortlock et al 



civ et ai.i 
J |201lJ). 



l200d : iPerez-Gonzalez et all I200S 
Similarly, the peak of cosmic star formation activity has 
been observed to occur at a redshift of z > 1 (see e.g ., 
iMadau et alll 19961 : ICimatti et aT]|2004l : ljuneau et aLll2005l) , 
and additionally seems to depend on galaxy mass: the most 
massive galaxies are among the first galaxie s to form (e.g., 



Thomas et"aH 120051 : IClemens et al. I 120061 : ICimatti et ail 
20061 ) and are roughly in place by a redshift of z ~ 1, where 



their number densities reach similar values as in the local 
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Universe (e.g.. [Conselice et al. I [20071 ). These early massive 
galaxies at z > 1.5 are found to have on average red rest- 
frame colours as expected for a maximally old stellar popu - 
lation at these very early epoc hs (|Griitzbauch et ai]|2011bT ), 
however, iBauer et al] i|2011br i show that 50 - 80% of these 
massive red galaxie s might harbour dusty star formation 
i|Bauer et al.ll2011bh . 

It is then clear that the main epoch of galaxy forma- 
tion occurs at z > 1, but it is still unclear which processes 
trigger or suppress star formation, and what are the most 
important factors in determining galaxy properties. Is it pri- 
marily galaxy mass (i.e. a galaxy's gravitational potential), 
is it galaxy merging on local scales or is it the large scale 
environment that drives galaxy evolution? 

Various processes connected to a galaxy's environment 
are expected to trigger or suppress star formation and in- 
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fluence galaxy assembly dGunn fc Got j 19721; Larson et al.l 



ll980l ; lMihos fc Hernquisdll994l ; lMoore et alj|l996h . The effl- 
ciency of these processes depends on the local galaxy density 
(e.g., for minor and major merging or harassment) or the 
density of the intra-cluster medium (e.g., for ram-pressure 
stripping). We might then expect to find a correlation be- 
tween local density and star formation rate, as it is in 
deed found in the local Universe (e.g., I Lewis et al.l 



Bal ogh et ai]|2004) and up to z ~ 1 (e . g..lElbaz et alj 



Cooper et all 120081 ; iPatel et all 120091 ; ISobral et all 1201 J ). 
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however with varying conclusions. At lower redshifts star 
formation rates of galaxies in clusters are found to decrease 
both with smaller cluster-centric dist ance and as a function 
of local projected galaxy density (e.g.. |Poggianti et al. 19991; 

1 Lewis et alj |2002| . and references therein). iBalogh et al.l 
(|2004 ) show that the correlation between local density and 
star formation rate extends beyond the dense cluster en- 
vironment to galaxies in the field. They argue that short- 
term local processes like galaxy merging at high redshift 
causes this dependence. At z ~ 1, some auth ors report a 
turn-around of the local SFR-density relation (|Elbaz et al.l 
120071 ; I Cooper et alj 120081 ) , corresponding to enhanced star 
formation activity in areas of high local density. Others (e.g., 
iPoggianti et al ]|200Sl ; IPatel et al.|[2009l l find that star forma- 
tion is suppressed in the highest density areas. This effect 
is observed in the cores of the most massive galaxy clus- 
ters a s early as z ~ 1.4 (|Lidman et al.l [20081 : iBauer et al.l 
l2011bh . while in less massive clusters at similar redshifts 
star formation in the cluster core was found to be enhanced 
|Havashi et~ai]|2010l; pTran et alj|20ich . 

Recentlv. lSobral et al.l (|201ll ) presented a study of SFRs 
over a wide range of environments at z ~ 1 and found 
that the median SFR and star-forming fraction increases 
with local density up to a critical surface density of E = 
10-30 Mpc -2 and decreases at higher densities, possibly 
partly reconciling the contradictory results of the above 
studies. However, lower SFRs and redder galaxy colours 
were also observed in general field sampl es outside the dense 
cores of galaxy clu sters up to z ~ 1.8 (|Quadri et al.l [201 ll : 
IChuter et al.ll20l"lT ). The controversy of the role of environ- 
mental processes in influencing the build-up of galaxy stellar 
mass at z > 1 still remains. 

So when and where does the relation between star 
formation and environment emerge? To witness the trans- 
formation from heavily star-forming to passively evolving 
galaxies we have to move to the crucial redshift domain at 

2 ~ 1.5 and above. In this epoch most of the stellar mass 
is built up and the formation of massive galaxies is stil l on- 
going (e.g.. iFontana et al.l [200^ : iMortlock et al.ll201ll ). To 
investigate the influence of galaxy environment on this pro- 
cess, a very deep, high spatial resolution survey, that allows 
for the detection of faint satellites, is crucial. In this study 
we utilise data from the GOODS NICMOS Survey (GNS), 
an extremely deep, near infrared Hubble Space Telescope 
survey, reaching unprecedented stellar mass completeness of 
log M, = 10 9,5 M ( 7) up t o z = 3. Using the same sample, 
iGriitzbauch et alj (|2011bT ) found that the colour-density re- 
lation is possibly reversed at z ~ 1.5 — 2 and is not detectable 
at z > 2. In the present study we investigate the behaviour of 
star formation rates as a function of local density as well as 
the possible differences between the SFR-density and colour- 
density relation introduced by the presence of strong dust 



extinction. Using a stellar mass limited sample also allows 
us to investigate the role of stellar mass relative to that of lo- 
cal environment without being biased towards star-forming 
galaxies, which is often the case for optically selected sam- 
ples. 

We present this data and the measurements of stellar 
masses, star formation rates and local densities in Section[2] 
The results are shown in Section [3] while we discuss and 
summarise our findings in Section 2] Throughout the paper 
we assume the standard ACDM cosmology, a flat universe 
with Qa = 0.73, Q,m = 0.27 and a Hubble constant of -Ho = 
72 km s" 1 Mpc -1 . 



2 DATA AND ANALYSIS 

In this section we describe the survey we use in this study, 
the GOODS NICMOS Survey (GNS), as well as the mea- 
surements of photometric redshifts, stellar masses, rest- 
frame colours, local densities and star formation rates. 



2.1 The GNS data 

The data used in this study is obtained with the GOODS 
NICMOS Survey (GNS). The GNS is a 180 orbit Hubble 
Space Telescope survey consisting of 60 single pointings with 
the NICMOS-3 near-infrared camera, with an imaging depth 
of three orbits per pointing. 

The pointings were optimised to contain the maximum 
number of massive galaxies (M» > 10 11 Mq) in the redshift 
range 1.7 < z < 3, identified in the two GOODS fie lds by 
their optical-to-infrared colours (|Conselice et alJ lboil). The 
survey covers a total area of about 45 arcmin 2 with a spatial 
resolution of ~0.1 arcsec/pixel, corresponding to ~0.9 kpc 
at the redshift range we use here (1.5 < z < 3). The target 
selection, su rvey characterist i cs an d data reduction are fully 
described in lConselice et aL j 201ll). For first s c ience results 
from the survey see e.g., iBuitrago et alj (|2008l).lBluck et all 
|2009l ). [Griitzbauch et al l (|2011rJ ) and lBauer et al.1 (|2011bh 

The GNS has a 5a limiting magnitude of Hab = 
26.8, which is significantly deeper than ground based near- 
infrared imaging of the GOODS fields done with ISAAC on 
the V LT, reaching a 5<r depth of Hab = 24.5 (|Retzlaff et al.l 
l2010h . Sources were extracted from the NICMOS H-band 
image and matched to the optical HST-ACS bands B, V , i 
and z, which is available down to a limiting magnitude of 
Bab = 28.2. The matching is done within a radius of 2", 
however the mean separation between optical and H-band 
coordinates is much better with ~ 0.28 ± 0.4", roughly cor- 
respon ding to the NICMOS resolution (see also lBauer et alj 
l2011bh . 

The resulting ff-band selected photometric catalogue 
covering the bands BVizH comprises about 8300 galaxies 
and is used to compute photometric redshifts, rest-frame 
colours and ste llar masses as describ ed in the following sec- 
tions (see also IConselice et al.ll201ll . for more details). 



2.2 Photometric Redshifts 

Photometric redshifts were obtained by fitting template 
spectra to the BVizH photometric data points using the 
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HYPERZ code dBolzonella et al.ll200(j). The method is de- 
scribed in more detail in iGriitzbauch et al.l l|2011bh . The 
synthetic spectra used by HYPERZ are c onstructed with the 
Bruz ual & Chariot evolutionary code |Bruzual fc Chariot! 
Il993h representing roughly the different morphological types 
of galaxies found in the local universe. We use five template 
spectra corresponding to the spectral types of E, Sa, Sc and 
Im as well as a single bur s t scen ario. The reddening law is 
taken from ICalzetti et al.l (|2000h . HYPERZ computes the 
most likely redshift solution in the parameter space of age, 
metallicity and reddening. The best fit redshift and corre- 
sponding probability are then output together with the best 
fit parameters of spectral type, age, metallicity, Ay and sec- 
ondary solutions. 

To assess the reliability of our photometric redshifts 
we compare them to available spectroscopic redshifts i n 
the GOODS fields (|Barger et al.ll2008l ; IWuvts et al.ll200Sl) . 
We matched the two catalogues to our photometric cata- 
logue with a matching radius of 2", obtaining 906 secure 
spectroscopic redshifts. The reliability of photometric red- 
shift measures is usually defined as Az/(1 + z) = (z spec — 
z p hot)/{l + z spec ). In the following we compare the median 
offset from the one-to-one relation between photometric and 
spectroscopic redshifts, {Az/(1 + z)), and the RMS scat- 
ter around this relation, o"a z /(i+z)- For the whole sample 
over the full redshift range we obtain a median offset of 
(Az/(l + z)) = 0.011 and a scatter of a Az/(1+z) = 0.061. We 
then investigate the performance of HYPERZ at different 
redshifts, at low redshift (z < 1.5) and in the redshift range 
of 1.5 < z < 3, which is the redshift range of the galaxy 
sample we use in this study. For the high redshift sample 
we obtain an average offset (Az/(1 + z)) = 0.06 and a RMS 
of ctaz/(i+2) = 0.10, with a fraction of catastrophic outliers 
of 20%, where catastrophic outliers are defined as galaxies 
with |Az/(l + z)| > 0.3, which corresponds to ~ 3 times the 
RMS scatter. Galaxies below z = 1.5 show a slightly lower, 
but still comparable scatter of (Jaz/{i+z) = 0.08, however 
the outlier fraction decreases dramatically to only ~ 2%. 
These are the values that are used in the Monte Carlo Sim- 
ulations described in Section 12.71 to account for the effect 
of the photometric redshift errors on the local density and 
SFR estimates. 



2.3 Stellar masses and rest frame colours 

The determinati on of stellar masses for our sample is 
fully d escribed in lConselice et all l|201ll) and lMortlock et al.l 
(2011). The stellar masses and rest-frame colours we use 
are measured by multicolour stellar population fitting tech- 
niques, based on the same catalogue used for the photo- 
metric redshift measurements. Spectroscopic redshifts are 
used if available. A large set of synthetic spectral energy 
distributions ( SEDs) is constructed from the stellar popula- 
tion models of iBruzual fc Charlotl <|2003l . BC03), assuming 
a Salpeter initial mass function. The star formation history 
is characterised by an exponentially declining model with 
various ages, metallicities and dust extinctions. The model 
SEDs are then fit to the observed photometric datapoints 
of each galaxy using a Bayesian approach. For each galaxy 
a likelihood distribution for the stellar mass, age, and abso- 
lute magnitude at all star formation histories is computed. 
The peak of the likelihood distribution is then adopted as 



the galaxy's stellar mass and absolute U and _B-band mag- 
nitudes, while the uncertainty of these values is given by the 
width of the distribution. We chose to compute rest-frame 
(U — B) colours, since the wavelength range of the U and B 
bands is covered best by the observed optical and H bands. 

While parameters such as age, e-folding time and metal- 
licity are not accurately fit due to various degeneracies, the 
stellar masses and colours are robust. From the width of 
the probability distribution we determine typical errors for 
our stellar masses, which are mainly due to uncertainties in 
the template fitting and photometric errors. There are ad- 
ditional uncertainties from the choice of the IMF, which are 
not taken into account here. We obtain a total random error 
of our stellar masses of ~0.3 dex, roughly a factor of two. 
A more detailed discussion of the stell ar mass uncertainties 
can be found in lConselice et all l|201lf) . 

We do not use any photometric data redder than the 
NICMOS //-band in our measurements of stellar masses. 
The reason for this is that essentially there is no data redder 
than the //—band which has the same fidelity and depth as 
the BVizH bands we use in this paper. The K-band data 
available from the ground is nowhere near as deep as the 
//-band NICMOS data. While we have IRAC data for our 
sources, we do not use these data due to the PSF issues and 
contamination from neighbouring galaxies. Furthermore, the 
rest-frame /3-band gives us a go od anchor for measur ing 
stellar masses, as is shown by e.g . .lBell fc de Jond l|200lf) at 
lower redshifts, and lBundv et all |2006l ) for higher redshifts. 

2.4 Completeness limits 

We compute the expected stellar mass completeness lim- 
its from the 5<r magnitude limit of the GNS (Hab = 26.8) 
and the mass-to- light ratios (M/L) of simple stellar popula- 
tions (SS P). For more de t ails o f this procedure we r e fer th e 
reader to lConselice et al.l (|201ll ) and lMortlock et al.l l|201ll ). 
The limiting stellar mass for a galaxy with a maximally old 
stellar population at z — 3 is M* = 10 9 ' 2 A/q. The influence 
of the presence of du st on the stellar ma ss completeness 
limit is investigated in iBauer et al.l (|2011bh . They find that 
a very high dust extinction of Ay > 3 magnitudes is neces- 
sary to remove low-mass heavily dust-obscured star-forming 
galaxies from our sample. It is therefore unlikely that our 
sample is significantly biased against dusty star-forming ob- 
jects. To make sure that our sample is not biased towards 
blue galaxies at high z we use a conservative stellar mass cut 
of log M 4 = 9.5 in the following. The colour dependent com- 
pleteness limits and mass functions of red a nd blue galaxies 
in the GNS are investigated in more detail by Mortlock ct al. 

Hon]). 

In this study we will focus on the redshift range of 1.5 < 
z < 3, which provides a co-moving survey volume of ~ 2.3 x 
10 5 Mpc 3 , minimising the effects of cosmic variance. The 
cosmic variance associated with a certain co-moving volume 
depends on the clustering strength of the objects and can be 
roughly estimated from the average galaxy number density 
and the expected variance of d ark matter halos at a certain 
redshift (|Somerville et al.ll2004f) . For the co-moving volume 
of the GNS within the redshift range of 1.5 < z < 3, and at 
the average stellar mass of galaxies in the survey, we expect 
the influence of cosmic variance to be less than 10%. The 
final galaxy catalogue we use in the following comprises 1289 



4 Griitzbauch et al. 



galaxies down to a stellar mass of log M» =9.5 within the 
redshift range of 1.5 < z < 3. 

2.5 Local densities 

The local densities w e use here are measured by 
iGriitzbauch et alJ <|2011bl ) using the same sample to perform 
a study of galaxy colours and blue fractions in different en- 
vironments. Two different approaches are considered in the 
above study, (1) the aperture density, based on galaxy counts 
in a fixed physical aperture and (2) the nearest neighbour 
density, based on the distance to the 3 rd , 5 th and 7 th near- 
est neighbour. Galaxies down to the completeness limit of 
log M» = 9.5 are taken into account here. For both methods 
a redshift interval of Az — ±0.25 is used to minimise the 
contami nation of our samp l e with f oreground or background 
objects. IGriitzbauch et all (|2011bl ) show that the 3 rd near- 
est neighbour density is best suited to identify extremes in 
the local density distribution, which is why we will use it as 
a local density estimator in this study. 

We give an overview of the method in the following. 
First, the distance to the third nearest galaxy, D3, within 
the redshift interval of Az — ±0.25 is computed for each 
galaxy. Secondly, we have to account for edge effects, which 
are a major issue due to the design of the GNS, since it 
does not have a continuous survey area. The area contain- 
ing the three nearest neighbours might not be fully covered 
for all galaxies, especially for galaxies close to the edges of 
isolated, non-overlapping pointings, since the field of view 
of a single pointing covers about 500 kpc at z > 1.5. To 
properly account for the loss of area due to survey edges 
we approximate the sampled area around each galaxy by 
the number of image pixels, which are actually covered by 
the observations, within the area 7rZJ§. Since the number of 
pixels within nDs (N p i xe i s (D3)) is directly proportional to 
the covered area, it can be used instead of ~kD\ to compute 
the surface density E3 = 3/N p i x< ,i s (Ds). This gives £3 in 
arbitrary units, which is ideal for our purposes, since we are 
only interested in relative densities and £3 is normalised by 
the median value in each redshift slice, (S3) a*. The nearest 
neighbour density in units of a relative over-density is then 
given by: 

(1 + S 3 ) = -^- (1) 

where £3 itself is the over-density. In the following we use 
log (1 + 5$) to distinguish between under-dense (log (1 ± 
S3) < 0) and over-dense (log (1 + S3) > 0) environments. 

2.6 Star formation rates 

The star formation rates (SFRs) were measured from the 
re st-frame ultr a violet (UV) luminosities, as fully described 
m iBauer et al.1 |2011bh . The rest-frame UV provides one of 
the most direct measurements of ongoing SFR, since the 
UV luminosity is directly related to the presence of a young 
and short-lived stellar population produced by recent star 
formation. However, UV light is very susceptible to dust 
extinction and a careful dust-correction has to be applied. 
The correction we use here is based on the rest-frame UV 
slope, i.e. a far-to-near UV colour. The method is described 



111 IBauer et all (|2011bh . presenting a detailed study the star- 
formation properties of the same sample we use here. We 
briefly describe the method in the following. 

We determine the SFRyv from the observed optical 
ACS/zgso-band flux density (with a 5 a limit of 27.5 in 
the AB system), which corresponds to the rest-frame UV 
luminosity around 2800A, spanning wavelengths of 2125 - 
3400A for z = 1.5 - 3 galaxies. First we derive absolute 
magn itudes using the kcorrect package (|Blanton fc Roweisl 
120071 . v4.2). The 2800A luminosity is then converted into a 
star formation ra t e assu ming a Salpeter IMF and using the 
iKennicutt et al. I i|l998T ) law: 

SFRuv (M© yr" 1 ) = 1.4 x 10" 28 Lasoo (ergs s" 1 Hz" 1 ) (2) 

A dust correction is determined for each galaxies uniquely 
by fitting model SEDs to the observed galaxy magnitudes 
from the optical t o the mid-infrared, fo l lowing the proce- 
dure described in IPerez-Gonzalez et al.l (|200Sf ) . The best- 
fitting template is used to obtain synthetic UV luminosities 
at 1600A and 2800A. The UV-slope /3 is then used to derive 
th e extinction at 2800A , A2800 following the law calibrated 
bv lCalzetti et al.l (|200Ch . 

The depth of the z-band data from HST-ACS allows 
us to reach limiting SFRs of 1.5 Mqyt' 1 at z = 1.5 and 
5 M©yr _1 at z = 3. These values are determined by using 
the 2-band limit plus a one magnitude dust correction to cal- 
culate the SFR across the redshift range used in this study. 
Of the 1289 galaxies within the redshift range of 1.5 < 2 < 3 
and stellar masses above M* = 1O 9 ' 5 M0, the majority of 
galaxies show signs of star formation, with only 1.5% not 
detected in the z-band. 

2.7 Monte Carlo simulations 

To estimate the reliability of the local density estimates we 
perform a se t of Monte Carlo s i mulati ons, using the same 
method as in IGriitzbauch et all (|2011bh . To asses the effect 
of the photometric redshift uncertainties on the local densi- 
ties and star formation rates we randomise the input redshift 
catalogue according to the Az/(1 + z) error obtained by the 
comparison with available secure spectroscopic redshifts in 
Section 12.21 We use the typical photometric redshift error 
Az/(1 + z) — 0.10 in our redshift range of 1.5 < z < 3 and 
assume a Gaussian distribution of errors, where the width 
of the distribution a corresponds to Az/(1 + z). For each 
galaxy a random value is selected within this distribution, 
which is then added to the measured photo-z. To account for 
scattering in and out of the redshift range we include the full 
GNS sample in the randomisation and compute the redshift 
uncertainty for each redshift range. We obtain an average 
Az/(1 + z) = 0.08 for galaxies with redshifts lower than our 
range of interest (z < 1.5) and an average Az/(1 ± z) — 0.1 
for galaxies with redshifts higher than the range we use in 
this study (z > 3). We deal with catastrophic outliers sep- 
arately by randomly adding much larger offsets for a per- 
centage of galaxies, corresponding to the catastrophic outlier 
fraction at the respective redshift (see Section [2. 2|l . Galax- 
ies with |Az/(l ± z)\ > 0.3 are treated as catastrophic out- 
liers here. The offsets are randomly picked from the interval 
0.3 < Az/(1 ± z) < 1, and added or subtracted from the 
original redshift. 

The randomisation procedure is repeated 100 times for 
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Figure 1. Dust corrected star formation rates (SFR, top) and specific star formation rates (SSFR, bottom) as a function of local density 
log (1 + (53). The small symbols are the observed data, while the big symbols represent the average of all Mont e-Carlo simulations. Each 
line is the result o f a single Monte-Carlo run. The red squares are lower redshift data from lLewis et ab (2002) (SFR, upper panel) and 
iPatel et all tolll ) (SSFR, lower panel), see text for details. 



each galaxy in the sample, resulting in a randomised photo- 
metric redshift catalogue, which is then used to recalculate 
the local densities and star formation rates for each of the 
100 Monte Carlo runs. From this we obtain 100 randomised 
local density estimates for each galaxy in the sample, from 
which we then compute the average local density uncertainty 
for the whole sample. The local density error for each indi- 
vidual galaxy is simply the standard deviation of all Monte 
Carlo runs. The individual errors are then averaged to ob- 
tain the average uncertainty of the local density estimator, 
which is A log (1 + <5 3 ) = 0.24. The SFR of each galaxy is 
recomputed using the change in luminosity distance caused 
by the change in photometric redshift in each simulation to 
obtain a new UV luminosity L2800, which is then used in 
equation [2] as above. 

The results of the Monte Carlo simulations are used 
throughout the paper to discuss the relations between SFRs, 
local density and stellar mass in the following sections, 
and are plotted in the respective figures. This allows us to 
demonstrate the uncertainties introduced by the photomet- 
ric redshift error and to asses the reliability of the observed 
trends. In all figures the original data points are plotted as 
small symbols, whereas each Monte Carlo run is represented 
by a solid line showing the average in bins of local density 
(figures [T] and [3]) or stellar mass (figure [4]) or colour (fig- 



ure 0, as further described below. The average of all runs 
in each bin is shown in the figures as big solid symbols with 
corresponding errorbars. We use the bi-weight estimator for 
locatio n and scale, following the definition of iBeers et al. I 
(1990), to obtain a robust estimate of the average and scat- 
ter in each bin. This method was originally devised for mea- 
suring the velocity dispersion of galaxy clusters, but can be 
generally applied to measure a robust mean and dispersion 
in a population following a non-Gaussian distribution with 
significant fractions of outliers as well as small sample sizes. 



3 RESULTS 

In the following we investigate the relations between star 
formation rates (SFR), specific SFR (SSFR) and local den- 
sity. The specific star formation rate is defined as the SFR 
per unit of stellar mass, M», i.e., SSFR = SFR / M*. It 
gives the relative importance of star formation with respect 
to the already existing stellar mass of a galaxy. In the fol- 
lowing sections we will first present the SFR-density rela- 
tion in three redshift bins between z = 1.5 and z — 3 over 
the whole stellar mass complete sample down to log M* = 
9.5 (Section 13. ip . We then investigate the connection to 
the colour-density relation found in the same sample by 



6 Griitzbauch et al. 



iGriitzbauch et all (|2011bl ) and the effect of dust extinction 
(Section 13. 2|) as well as the role of stellar mass in the SFR- 
density relation (Section 13. 3[) . For this purpose we split the 
sample in low and high quartile of stellar mass and local den- 
sity respectively. Finally, in Section [3T4l we show the influence 
of the very local environment (number of neighbours within 
30 kpc) on star formation rates and rest-frame colours. 



3.1 The SFR-density relation 

Figured] shows SFR (top) and SSFR (bottom) of all individ- 
ual galaxies in our sample down to the completeness limit 
of log M„ — 9.5 as a function of relative local density. The 
observed data points are shown as small symbols. No strong 
trend between SFR or SSFR and local density is present in 
our data. 

To asses the significance of the result we also plot the 
results of the Monte Carlo simulations. Each solid line in 
Figure [T] represents one Monte Carlo run, averaged in bins 
of local density of A log (I + 83) = 0.6. The average SFR and 
SSFR of each bin is plotted against the average log (1 + 53) 
for each run. The average and scatter over all runs is shown 
as big magenta circles (SFR) and blue triangles (SSFR) in 
the respective panel. The three columns show the results in 
three redshift bins. To show possible evolution in the above 
relations, we have divided our redshift range of 1.5 < z < 3 
into three bins with a bin size of Az = 0.5. The first column 
then corresponds to the redshift range of 1.5 < z < 2, the 
second to 2 < z < 2.5 and the last to 2.5 < z < 3. 

The top panel shows the SFR-density relation in the 
three redshift bins. We do not see a strong general trend 
for increasing or decreasing SFR with local density at any 
redshift, apart from a lower average SFR in the highest 
local density bin at the lowest redshift (1.5 < z < 2). 
The difference in average SFR between over-densities and 
average/under-dense areas is about a factor of three. A K-S 
test shows that it has a significance of just over 2a. The av- 
erage specific star formation rates (SSFRs) show a similar 
behaviour with local density as the average SFRs. There is 
a decrease of SSFR with local density in the lowest redshift 
bin (1.5 < z < 2) with a significance of ~ 2a. 

For comparison re asons we hav e plott ed the local SFR- 
density relation from iLewis et al.l (l2002Tl and the SSFR - 
density relation at z ~ 0.6 — 0.9 from Patel et al.l l|2011T l. 
Both data-sets are converted into relative over-densities by 
normalizing the absolute densities by the median density of 
the whole sample quoted in each study. We have then scaled 
both, the median SFRs and the median SSFRs by a factor 
of 10 to keep a convenient scale in the figure. Both data-sets 
are plotted as red squares in the panel of the first redshift 
bin in Figure [1] The slope of both, the SFR- and SSFR- 
density relation is considerable steeper in the lower redshift 
samples. This will be further discussed in section [4] 

We do not find any correlation between SFR or SSFR 
and local density at z > 2. This might partly be due to 
the increasing scatter between the single Monte Carlo sim- 
ulations. Due to the scatter at low and high densities we 
cannot rule out the presence of a change in average SFR 
of a factor of up to three (at the 2a level) between the ex- 
tremes of local density. Larger samples will be needed to 
address this in more detail. 



3.2 The correlation to restframe (U — B) colour 
and the effect of dust extinction 

Galaxy colour is often used as a proxy for star-formation 
activity and history, however heavy dust extinction can con- 
ceal substantial amounts of current star formation leading 
to red res t-frame colours indi cative of an old stellar popula- 
tion fe.g.. lBauer et alll2011bh . Our data allows us to assess 
if the extrapolation from colours to star-formation activity 
is valid at z > 1.5 and how the two properties correlate with 
local environment. 

So how does the SFR-density relation compare 
to the colour-density relation in the same data-set? 
IGriitzbauch et ah! (|2011bh find that overall there is no strong 
correlation between (U — B) colour and local density at 
z > 1.5, however, the highest over-densities tend to be in- 
habited by a higher fraction of blue galaxies than interme- 
diate and low density environments. If red colour implies 
less star formation activity, then this trend seems to be 
opposite to what we find here for the SFRs and SSFRs, 
which both tend to be lower in the highest over-densities, 
however only in the lowest redshift range we probe here 
(1.5 < z < 2). This apparently contradictory behaviour 
of colours and SFRs in the highest density bin can be ex- 
plained by the positive correlation between SFR and colour 
at z > 1.5: the redder the galaxies the higher their av- 
erage SFR. This trend is illustrated in Figure [2] which 
shows the average SFR and SSFR as a function of (U — B) 
colour. The SFRs of red galaxies are boosted by the dust- 
correction which strongly correlates with rest-fr ame (U — B) 
colour. This trend is discussed in more detail by I Bauer et al.l 
|2011bl ). They show that the extinction values increase from 
A v ~ 1 at (U - B) ~ 0.5 to Ay > 3 at (U — B) > 1. This 
dependence is supported by both the UV-slope /3 and the 
24^im-detection fraction, which both increase with (U — B) 
colour. 

The specific star formation rates decrease with colour, 
as expected, since both colour and SSFR trace the relative 
importance of the current star formation with respect to 
the total stellar population of a galaxy. This means that 
the high SFRs of red galaxies are due to the on average 
hi gher stellar masses of red g alaxies (see also the discussion 
in IGriitzbauch et al.1 (fcullbl )'). 

The dust-correction is a major source of uncertainty 
in the determi n ation o f SFRs from UV light, as discussed 
in lBauer et al.l (|2011bh . The dust-correction might bias the 
SFRs of red galaxies towards higher values, in case their 
red UV colour is due to an evolved stellar population rather 
than dusty star formation. This possible bias would nat- 
urally introduce a positive correlation between colour and 
SFR in the sample, such that red galaxies have higher SFRs. 
Any existing trend between colour and local density would 
then result in a similar relation between SFR and local den- 
sity. Note, however, that the trend of SFR and SSFR de- 
creasing with local density is present in both dust-corrected 
and uncorrected SFRs at 1.5 < z < 2. As stated above, al- 
though the dust correction derived from the UV-slope has 
a large uncertainty, it agrees well with the 24pm-detection 
fraction, a complete ly independent indicator of dust-content 
l|Bauer et al.ll2011bl ). 
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Figure 2. Star formation rate (SFR, top) and specific star formation rate (SSFR, bottom) as a function of rest-frame colour (U — B). 
Same symbols as in Figure [l] 



3.3 The role of stellar mass in the SFR-density 
relation 

It has been suggested by various authors that the de- 
pendence of galaxy properties on local density varies 
strongly with galaxy stellar mass (Cassata c t al.l 2007I; 
Bolzonella et all |2009 |; iTasca et all 120091: iPannella et all 



20091; 



Iovino et alj l20ld ; 
20111 ). 
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Sobral et all 120111 ). To investigate the impact of stel- 
lar mass on the SFR-density relation in our sample we now 
look at the SFR-density relation at "fixed" stellar mass. For 
comparison reasons we also investigate the relation between 
SFR and stellar mass at fixed local density. For this purpose 
we split the sample in high and low quartiles of stellar mass 
as well as high and low quartiles of local density. Note that 
due to our limited sample size we show the results for the 
whole redshift range 1.5 < z < 3 in these figures. To avoid 
confusion we plot only the Monte-Carlo-simulated data and 
average and dispersion in each bin of local density or stellar 
mass, as described in the beginning of Section [3] 

First we investigate the behaviour of average SFR and 
SSFR with density in the high and low quartile of the stellar 
mass distribution. The sample is split in the low quartile at 
log M» < 9.7 and in the high quartile at log M* > 10.5, 
containing by definition 25% of the sample respectively. Fig- 
ure [3] shows the SFR-density (top panel) and SSFR-density 
(bottom panel) relation in high (red, solid triangles) and 



low (blue, solid squares) quartile of stellar mass. We find no 
trend of SFR or SSFR with density within the errors, either 
for low mass or high mass galaxies. 

Figure U shows the average SFR and SSFR as a func- 
tion of stellar mass. The colour coding is the same as before, 
but now the sample is split in high and low quartiles of lo- 
cal density, including galaxies at log (1 + 5s) > 0.25 and 
log (I + 53) < 0.25 in the high and low quartile respectively. 
There is a strong correlation between SFR and stellar mass 
at all densities. The average SFR increases monotonically 
with stellar mass until about log M* ~ 11, where the re- 
lation flattens out at an average value of (SFR) ~ 100 Mq 
yr~ . The same can be seen in the relation between SSFR 
and stellar mass in the bottom panel of Figure [4] There is a 
strong anti-correl ation between SSFR and stellar mass (as 
also discussed by iBauer et al.l (|2011bh ). analogous to the 
strong correlation betwee n colour and stella r mass in the 
same sample discussed in iGriitzbauch et al l (|2011bl ). The 
correlation steepens towards the high mass end, as a result 
of the flattening of the SFR-stellar mass relation discussed 
above. This suggests that the relative importance of star 
formation in the most massive galaxies is already declining 
at very early times [z ~ 3), possibly connected to AGN- 
feedback. 

To investigate this possibility we compare the SFRs of 
likely AGN candidates in our sample selected via X-ray and 
IR-excess by Weinzirl (2011, submitted) with galaxies of 
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Figure 3. Monte Carlo simulation of SFRs (top) and SSFRs 
(bottom) as a function of local density log (1 + £3) in the high 
(red, triangles) and low (blue, squares) quartile of stellar mass. 



similar stellar mass that show no signs of hosting an AGN. 
We do not see any trend for a different distribution of SFRs 
in AGN hosts at fixed stellar mass, neither is the flatten- 
ing of the SFR-mass relation affected by the presence of 
AGNs in our sample. This is in agreement with a study 
of X-ray selected A GN host galaxies at 0.4 < z < 1.4 by 
iBundv et al.l ((2008), who find that in the reddest and most 
massive host galaxies, AGN may not be directly responsi- 
ble for quenching star formation. The flattening could then 
simply be caused by the quick exhaustion of the galaxies' 
cold gas supply following a period of an extreme star-burst 
in the most massive galaxies at even higher redshifts that we 
probe here (z > 3), p ossibly connected to galaxy merging 
i|Conselice et aJj|2008f ). 

The large scatter between the individual simulations at 



^ log (l + « 3 ) > 0.25 
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log M* 



^ log (1+6,) > 0.25 
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Figure 4. Monte Carlo simulation of SFRs (top) and SSFRs 
(bottom) as a function of stellar mass log M* in the high (red, 
triangles) and low (blue, squares) quartile of local density. 



the low and especially the high density end is visible in all 
figures. This scatter is caused by the low number of objects 
in extreme densities. The scatter is also much larger at high 
stellar masses, due to the intrinsic scarcity of the most mas- 
sive galaxies. The large scatter at high stellar masses above 
log M* > 10.5 could hide a difference in SFRs between 
the density extremes of a factor of ~ 5 at the 2a level. In 
the low mass quartile the dispersion is much smaller, but 
still possibly hiding a factor of ~ 2.5. We conclude that the 
present sample is not well suited to investigate the behaviour 
of the SFR-density relation at different stellar masses due 
to the limited cosmological volume of the survey, which 
is especially problematic for high mass galaxies. However, 
we can exclude the presence of a strong SFR-density re- 
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satellites within 30 kpc 

Figure 5. SFR and SSFR as a function of number of close neighbours. The original data is shown as small symbols, whereas the results 
of the Monte Carlo simulations is overplotted as lines (one for each individual simulation) and big symbols (overall average). 



lation comparable to the local relation for galaxies below 
log Af* = 10.5. 

3.4 The role of the very local environment within 
30 kpc 

Which mechanism might be mainly responsible for the on 
average lower SFRs at the highest over-densities observed in 
our sample at z < 2? One possible origin would be the ear- 
lier and faster formation o f more massive galax ies in higher 
density environments fe.g.. lThomas et al.ll2005t ). This would 
imply that there should be a similar relation between lo- 
cal density and stellar mass, such that high-mass galaxies 
tend to be located in high over-densities. Another possible 
mechanism could be galaxy merging, which can trigger a 
short-lived burst of star formation exhausting a galaxy's gas 
reservoir and subsequently suppressing star formation, how- 
ever, th e time-scales in v olved in this process are very un- 
certain. IConselice et al.l l|2008l ) have shown that the galaxy 
merger rate peaks at z ~ 2, using galaxy morphologies 
and asymmetries in th e light distribu t ion (C AS parameters, 
IConselice et~aH (|2003h V iBluck et al.l ((2009) using the mas- 
sive galaxies from our sample (log M, > 11) have found 
that close neighbour counts within a radius of 30 kpc are a 
good tracer of the galaxy merger rate up to z ~ 3, leading to 
ve ry similar result s as th e morphology based merger rates 
in IConselice et alj (2009). Furthermore they find that the 



peak in merger fraction happens earlier for massive galaxies 
compared to less massive ones. This enhanced merging ac- 
tivity could then induce a higher star-formation activity in 
massive galaxies dependent on redshift and on the number 
of close neighbours. 

To test the above hypotheses we now investigate the 
behaviour of SFR and SSFR, stellar mass and rest-frame 
colour as a function of the number of neighbours within a 
radius of 30 kpc. For this purpose we count the number 
of galaxies located within a 30 kpc co-moving radius and 
a certain redshift interval for each galaxy in our sample. 
To minimise foreground and background contamination the 
same redshift interval of Az = ±0.25 as in the local density 
computation is used here. 

Figure[5]shows the SFR and SSFR as a function of num- 
ber of neighbours within 30 kpc. The original data are plot- 
ted as small symbols (discrete values), whereas the Monte 
Carlo simulated data is shown as solid lines for each individ- 
ual randomisation and as big symbols for the average and 
spread of all runs. Interestingly, we see the same trend of 
lower average SFR at higher local densities as for log (1+^3), 
but here the SFR is decreasing more gradually with number 
of neighbours. The trend is most visible in the lower red- 
shift bin and weakens with increasing redshift. Notice that 
there are only very few objects with three or more neigh- 
bours. We estimate the significance of this trend through 
a K-S test between the SFRs and SSFRs of galaxies with 
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Figure 6. Stellar mass and rest-frame colour as a function of number of close neighbours. The original data is shown as small symbols, 
whereas the results of the Monte Carlo simulations is overplotted as lines (one for each individual simulation) and big symbols (overall 
average) . 



two or more neighbours and those of galaxies with less than 
two neighbours. This suggests that the difference in aver- 
age SFR is significant at the 3<r level in the lowest redshift 
bin (1.5 < z < 2) and at up to z ~ 2.5. The average 
SSFR of galaxies with two or more neighbours are signifi- 
cantly lower (3a) than the SSFRs of galaxies with less than 
two neighbours up to z < 2.5. 

Is this trend due to a variation of stellar mass or colour? 
Figure [6] shows the dependence of stellar mass and rest- 
frame colour on the number of close neighbours. We see no 
significant dependence of average stellar mass or (U — B) on 
the number of close neighbours. The lower average SFR for 
galaxies with close neighbours then seems to be not a stellar 
mass induced trend but is genuinely connected to a galaxy's 
environment. 



4 DISCUSSION 

The major conclusion we draw from this study is the pres- 
ence of a strong dependence of SFR and specific SFR on 
galax y stellar mass, whic h is already detectable at z ~ 3 (see 
also lfBauer et al.ll2011bl ). The average SFR increases with 
stellar mass reaching values of SFR ~ 100 at log M„ ~ 11. 
We find a flattening of the SFR-stellar mass relation at the 
highest stellar mass end (log M» > 11), possibly connected 



to AGN feedback. However, we do not see significantly dif- 
ferent SFRs or specific SFRs in AGN candidates at fixed 
stellar mass. The average specific SFR decreases with stel- 
lar mass, indicating that star formation activity in massive 
galaxies is already beginning to decline at z ~ 3. 

In contrast to the strong correlation between SFR and 
stellar mass, overall we do not see a strong correlation be- 
tween SFR a nd local density, a s it is seen in the local 
universe (e . g.. lEewis et al.l 120021 ) and up to z ~ 1 (e.g., 
iPatel et all 120091 ). However, we find evidence for a possi- 
ble decrease in average SFR, as well as SFR per unit mass 
(SSFR), at the highest local densities at 1.5 < z < 2. Galax- 
ies in over-densities of more than a factor of five relative to 
the average local density exhibit on average lower SFRs as 
well as lower SSFRs than galaxies in average and under- 
dense regions. The difference amounts to a factor of 2-3 
which is less than the diff erence in the SFR-density rela- 
tion in the local Univers e llLewis et al] 120021: iBaloeh et aD 
l2004h and up to z ~ 1 (|Patel et all |2009t l201lh . who find 
a difference of at least a factor of five over a density range 
of two orders of magnitude. We do not detect any differ- 
ences in the SFR-density relation if the sample is split in 
low and high quartiles of stellar mass, while the SFR-stellar 
mass relation is equally strong at both high and low local 
densities. 

Our data is also compatible with no correlation between 
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SFR and local density, as was also foun d in a study of Lyman 
Break Galaxies (LBGs) at z ~ 3 by iBouche &: Lowenthall 
(2005). Note, however, that a sample of LBGs naturally in- 
cludes only blue, star-forming galaxies. As the authors point 
out, such a sample might not allow for the detection of a 
SFR-density relation present in a mass-complete sample. 
Similarly, the SFRs for star- forming galaxies at z ~ do 
not v ary strongly with local density (see e.g., iBalogh et al.l 
|2004) . 

Compared to the intermediate redshift range of z ~ 1 — 2 
our findings are not consistent with the in creasing SFRs with 
den sity in field ga l axies reported by e.g-. lElbaz et al.l (|2007T ) 
and lCooper et all ((2008) . Higher star- formation in the most 
dense environments w as also found in high redshift galaxy 
clusters at z = 1.46 l|Havashi et al.l l201Ch and z = 1.62 
l|Tran et alj|2010h . Other studies of a massive galaxy cluster 
at z ~ 1.4 ho wever show the opp osite trend of quenched 
star for mation (iBauer et aLlfeoilal ) and older stellar popu- 
lations |Lidman et al.ll2008l ) of galaxies in the cluster centre . 
Simi lar results we r e foun d recently by IChuter et all (|201lh 
and iQuadri et al.l (|201ll ) studying the colour-density and 
SFR-density relations at z ~ 1 — 2 in more general, field 
dominated samples. Both studies use the UKIDSS (UKIRT 
Infrared Deep Sky Survey) Ultra Deep Survey (UDS) and 
slightly different but similar techniques of measuring local 
densities (fixed aperture and nearest neighbour densities). 
The above authors show that galaxies in higher local den - 
sities have on average redder c olours JChuter et alj l201ll ) 
and lower star-formation rates (|Quadri et al.ll201ll ) out to 
2 ~ 1.8, even at fixed stellar mass. This is consistent with 
the redshift range in which we find on average slightly lower 
SFRs in the most over-dense environments in this study. 

One issue in comparing galaxy properties with local 
environment might be the different ranges of absolute lo- 
cal densities probe d by different studies, as pointed out by 
ISobral et all fcoill) . They use narrow band Ha imaging over 
a wide area to show that the median SFRs increase with lo- 
cal densities in environments that correspond to the field 
and into cluster outskirts, but drop sharply towards cluster 
core densities. Depending on the range of absolute densities 
and total mass of the structures that are probed by a given 
study, the result could vary from increasing over constant 
to decreasing SFRs with local density. To estimate the large 
scale density or total mass range of structures we are likely to 
probe in our s ample we app l y the cluster abundance func- 
tion found by iBahcall et all (|2003l ) to our survey volume. 
The abundances of structures of different mass is most likely 
evolving with cosmic time, however, it is not well constrained 
from observations at hig h redshift. We there fore use the lo- 
cal cluster abundances of lBahcall et al.l (|2003l ) to get a rough 
estimate. We find that we can expect only a few structures 
of a total mass within 0.6h -1 Mpc of M . 6 > 1O 14 M to be 
present in our sample. This suggests that in this study we 
are confined to the field and group environment and that 
it is very unlikely that we are probing environments with 
the highest absolute densities like massi ve galaxy clu s ters in 
this study. According to the results of ISobral et al.l (|201ll ) 
we should then find on average increasing SFRs with local 
density, which we do not see in our sample. 

We also investigate the influence of the presence of close 
neighbours within a radius of 30 kpc on star formation rates 
and find that galaxies with more close neighbours show lower 



average SFRs as well as lower SSFRs up to z < 2.5. This is 
consistent with what is found for the 3 rd nearest neighbour 
densities above, but the relation between SFR and num- 
ber of close neighbours has a higher significance than the 
SFR-density relation. We also see a more gradual decline of 
the average SFR with number of neighbours than with lo- 
cal density. The relation between SFR and number of close 
neighbours is not caused by variations in average stellar- 
mass. Star formation quenching induced by the very local 
environment in the highest density regions, possible due to 
galaxy merging, seems to be present in our sample at least 
up to z ~ 2 and possibly up to z ~ 2.5. 

Another possible explanation for our results is that lo- 
cal density traces the dark matter halo mass of the large 
scale structure in which a galaxy resides. The halo mass is 
expected to determine the time-scale and amount of cold 
gas accreted by a galax y, as predict e d by the hydrodynam- 
ical simulations of e.g.. iKeres et al.l (|2005l ). who model the 
evolution of cold gas accretion and star formation rates over 
cosmic time. Since the cold accretion fraction depends on 
local density - due to typically higher halo masses in high 
density areas - they also predict decreasing SFRs with lo- 
cal density. The typical density at which the SFR starts 
to decrease progressively shifts to higher densities at higher 
redshift and is detectable up to z ~ 2, but not at z ~ 3. 
This prediction is consistent with the low average SFR at 
the highest densities we see at z ~ 1.5 — 2. 



5 CONCLUSIONS 

We study the relationship between star formation, local den- 
sity and stellar mass at a redshift of 1.5 < z < 3 using data 
from a very deep, near-infrared HST-survey, the GOODS 
NICMOS Survey (GNS). We find the following three main 
conclusions: 

1. Star formation rates strongly depend on galaxy stellar 
mass over the whole redshift range of 1.5 < z < 3 we study 
here. The correlation between SFR as well as specific SFR 
and stellar mass does not depend strongly on environment 
and is similar at all local densities. 

2. The influence of the local environment on star forma- 
tion rates at z > 1.5 is much weaker than in the local uni- 
verse (and up to z ~ 1). The average SFR does not decrease 
monotonically with density, but is suppressed only in the 
highest density regions, and only up to z ~ 2, where over- 
densities of a factor of > 5 have on average lower SFRs (by 
a factor of ~ 3) than average and under-dense areas. This 
is in good agree ment with expect ations from hydrodynami- 
cal simulations (|Keres et al.ll2005T ). which predict that lower 
SFRs at the highest over-densities should be detectable up 
to z ~ 2. 

3. The presence of close neighbours within 30 kpc has a 
similar effect on the average SFRs as high 3 rd nearest neigh- 
bour densities. Overall both local density indicators show 
the same trend, but the decrease of average SFR in galaxies 
with several close companions is more significant than the 
decrease of average SFR at high relative over-densities. This 
suggests that the very local environment has an influence on 
the star formation activity in early galaxies possibly through 
the mechanism of galaxy merging. 
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In conclusion our data supports the emerging picture where 
stellar mass largely drives galaxy evolution. Local environ- 
ment seems to have little impact on the star formation rates 
of galaxies at z > 1.5, apart from the most dense environ- 
ments, possibly the most massive structures, which already 
had time to build up and exert their influence on galaxies 
after a lifetime of the Universe of only ~3 Gyrs. 
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